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This letter describes a novel derivation of general relativity by considering the (non)self- 
consistency of theories whose Hamiltonians are constraints. The constraints, from Hamilton's equa- 
tions, generate the evolution, while the evolution, in turn, must preserve the constraints. This 
closure requirement can be used as a selection mechanism for general relativity starting from a 
very simple set of assumptions. The configuration space is chosen to be a family of 3 x 3 positive 
definite symmetric matrices on some bare 3-manifold. A general Hamiltonian is constructed on 
this space of matrices which consists of a single constraint per space point. It is assumed that this 
constraint looks like an energy balance relationship. It will be the sum of a 'kinetic' term which is 
quadratic and undifferentiated in the momenta, and a 'potential' term, which is any function of the 
configuration variables. Further, the constraint must be a scalar under the linear group, the natural 
symmetry group of the configuration space. This inexorably leads to the ADM Hamiltonian for 
general relativity. Both the space of Riemannian geometries (Wheeler's superspace), and spacetime 
are emergent quantities in this analysis. 
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One can write general relativity in Hamiltonian form 
by choosing the configuration space as Wheeler's super- 
space, the space of 3-dimensional Riemannian metrics, 
gij, modulo diffeomorphisms. The Hamiltonian is Q 
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where ix 1 ' 3 , a symmetric tensor density, is the momen- 
tum conjugate to g^, R is the scalar curvature of g^, g 
is the determinant of <fc, and N and iVj, the lapse and 
shift, are Lagrange multipliers multiplying the Hamilto- 
nian and momentum constraints 
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The evolution equations preserve the constraints. One 
gets 
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where C is the Lie derivative in the direction of the shift. 

In 1976, Hojman, Kuchaf, and Teitelboim, 0], showed 
that if one considered a constrained Hamiltonian on 
superspace, whose Poisson bracket algebra mimicked 
Eas. l|4l5(l . then one was led to the ADM Hamiltonian, 
Eq.Q. They imposed these specific closure requirements 



because they wished to guarantee that they would recover 
spacetime and 4-diffeomorphism invariance. 

In a recent article, |3|, we showed that this was redun- 
dant. We chose superspace as the configuration space, 
left the momentum constraint as it was, but considered 
a (fairly) general Hamiltonian constraint. We assumed 
that the constraint was the sum of a kinetic term which 
was quadratic in the undifferentiated momenta and a po- 
tential term which was some general function of the met- 
ric. The requirement that the constraint algebra closed, 
without insisting on any particular form, was enough to 
recover the ADM Hamiltonian. 

One can do even better. In jj| I considered the ADM 
Hamiltonian without the momentum constraint. The 
Hamiltonian constraint is the only constraint and the 
evolution equation for it is just Eq.Q) without the Lie 
derivative term. The right hand side no longer weakly 
vanishes, we have a secondary constraint 
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Following Dirac, pj, we know we have to add this back 
to the Hamiltonian, with a Lagrange multiplier, and thus 
we recover the standard Hamiltonian. 

I wish to push this further. The obvious thing to do 
is to merge these two results and consider a single scalar 
constrained Hamiltonian on superspace. I intend here to 
go much further. I will abandon superspace, and replace 
it with the space of matrices with only the linear group as 
the symmetry. The requirement that the constraint al- 
gebra closes will force me back to 3-covariance, the ADM 
Hamitonian, and, eventually, to spacetime. 

I choose as my configuration space a bare 3- 
dimensional manifold on which I put a family of 3 x 3 
positive-definite symmetric matrices. I will write them 
as M a b- I can compute at each space point an inverse 
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matrix which I denote by M ab . These can be distin- 
guished because they behave differently under rigid ro- 
tations. The momenta are in the dual of the tangent 
space to the configuration space so I write them as P ab 
because they should transform just like the M ab, s, multi- 
plied by the determinant of the transformation. I intend 
to construct a scalar Hamiltonian. To do this I only need 
ensure that the 'up' indices match the 'down' indices. 

I further assume that the 'kinetic' part, K, of the 
Hamiltonian is quadratic in the undifferentiated mo- 
menta. There are only two ways we can sum the mo- 
menta. Therefore we have 

K = AM ab M cd P ac P bd - | [M ab P ab ] 2 , (7) 

where A and B are, as yet, arbitrary constants. I assume 
that the Hamiltonian also has a potential term Q which is 
some function on the space of matrices. My Hamiltonian, 
therefore, has the form 
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where N is a Lagrange multiplier. 

I will now consider the possible choices for the function 
Q. The simplest term that Q could have is an undiffer- 
entiated scalar, CM ab M ab . This works well, it allows the 
final theory to have a cosmological constant. I will put it 
aside for the time being and seek a more complicated ex- 
pression. The next simplest term is a term that is linear 
in the first derivative of M, M abyC . Such a term cannot be 
used because the upstairs indices come in pairs so cannot 
be used to sum over the three downstairs indices. There- 
fore the first nontrivial term in the potential must be at 
the next order. 

This term will consist of terms which are linear in the 
second derivatives of the matrices. It turns out that there 
are only two such 

Q = CM ab M cd M acM - DM ab M cd M abiCd (9) 

From a dimensional argument we should also include 
terms which are quadratic in the first derivative such as 



Q' = M ab M cd M ef M abiC M, 
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There are five such terms. Therefore the Hamiltonian 
to this order will have nine constants, two in the ki- 
netic term, seven in the potential term. The logic now 
is quite straightforward; one works out Hamilton's equa- 
tions, applies them to the constraint and sees whether 
there is any possibility that one gets no, or especially 
simple, secondary constraints for any particular choice of 
the constants so that the final constraint algebra closes. 

It turns out that one can simplify this problem hugely 
by linearizing everything. This means assuming that the 



matrices M ab can be written as M ab = I + em ab , where / 
is the identity matrix, e is some small parameter and m ab 
is of order unity. The momentum P ab — y / ep ab . I do not 
assume that N is small. In turn, this means that we can 
drop the terms quadratic in the first derivative and dis- 
card many of the terms arising in Hamilton's equations. 
Therefore I assume the following Hamiltonian 
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This has four constants, (^4, B, C, D). We can set two of 
them to, say, 1 by an overall scaling of the Hamiltonian 
and absorbing a constant into the definition of the mo- 
mentum. One should not do this immediately because 
the ones one picks may turn out to vanish. 

Now I work out the evolution equations (to leading 
order) and calculate the time derivative of the constraint. 
Hamilton's equations are: 
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plus higher order terms. Computing the evolution equa- 
tion for the constraint gives 
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I wish to consider the first term on the right hand 
side of Ea. ((T3|l . We have three choices: 1) A = 0: This is 
deeply boring, we end up with no dynamics, all the terms 
eventually vanish. 2) C = 0: This leads to a well-known 
and viable theory, sometimes called 'strong gravity' 0. 
In this theory, we have a nontrivial kinetic term in the 
Hamiltonian, but the potential term vanishes. In addi- 
tion, we have no need for the momentum constraint. 3) 
A 7^ 0, C 7^ 0: from rescaling, we can set A = C = 1. 
This is the most interesting possibility, it forces the ap- 
pearance of a secondary constraint 
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(the linearized momentum constraint!) on the theory. 

This must now be added back to the Hamiltonian, I 
write this in the form — 2N a P c j J b , where N a is the new 
Lagrange multiplier. The altered Hamilton's equations 
have to be derived and the time derivative of the con- 
straints calculated anew. The evolution equation for the 
new constraint is especially interesting. It reads 

J^Pf = (C - D)M a °M e fN, egf , (15) 

where we have dropped a term, +M e f Pf^Nf, which van- 
ishes when the constraints hold. The right hand side of 
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Eg. (|15fl must vanish (otherwise we get a neverending se- 
quence of constraints). Therefore we require C = D. 

The second term on the right hand side of Eq. (|13|) has 
still to be dealt with. If BC + 2AD - 3BD ^ then we 
get 



M ef ^N 2 [M ab P ab ] 



=>■ [M ab P ab ] = constant. 

(16) 

This leads to two different theories. One is GR in the 
constant mean curvature gauge, the other is a confor- 
mally invariant theory which we call conformal gravity 
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In this letter, I am much more interested in the other 
choice, i.e., where BC + 2AD - 2>BD = 0. Using A = 
C = D = 1, we immediately get B = 1. Hence, my 
selfconsistcnt Hamiltonian becomes 
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This is nothing more than the linearized ADM Hamilto- 
nian! 

Now I check whether the constraint algebra still closes 
if the higher order terms in Hamilton's equations are in- 
cluded. Not surprisingly, it does not work. Next, I add 
back in the five terms which are quadratic in the first 
derivatives with arbitrary coefficients and discover that 
there is a unique choice of coefficients which makes ev- 
erything close. I need to add 

A H = M ab M cd M ef (M de . c M abJ - M de , c M af>b 
- ~M ab>e M cdJ - ^M de , a M bf , c + jM de , a M c/) (J18) 
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and I get a fully selfconsistent constrained Hamiltonian 
on the space of matrices. No other combination works. 

Each individual term in the Hamiltonian is only invari- 
ant under the linear group. However, the total Hamilto- 
nian is invariant under a much larger group, the group of 
local 3-diffeomorphisms. This is because the combination 
of terms that makes up the potential term is scalar cur- 
vature, written in terms of first and second derivatives of 
the metric. Therefore I can promote the matrices to the 
space of positive-definite symmetric tensors. I can iden- 
tify the 3-matrices with a 3-metric and recognise that the 
configuration space is isomorphic to superspace. 

This theory has arbitrary Lagrange multipliers, N and 
N a . Therefore the evolution from a given set of ini- 
tial data, viewed as a curve in superspace, is under- 
determined. This is Wheeler's 'many fmgured time'. 
Each solution curve on superspace can be used to gen- 
erate a pseudo-Riemannian manifold by the rule g = 



— l/N 2 ,go a — N a . One finds that each 4-manifold con- 
structed in this fashion satisfies the Einstein equations 
and that the solution curves arising from a given set 
of initial data can be mapped into one another by 4- 
diffcomorphisms. Therefore, starting from a constrained 
Hamiltonian on the space of matrices, there is a straight 
path that leads to three dimensional Riemannian space, 
four dimensional spacetime, the Lorentz group, and the 
Einstein equations. 

One obvious question that may be asked is why this 
choice of a constrained Hamiltonian? The naive answer 
is that it works. A somewhat more sophisticated an- 
swer may be found by making a Legandre transforma- 
tion from the Hamiltonian to a Lagrangian. Further, a 
parametrised action can be constructed which will have 
the classic square root form of the Jacobi action (sec 8]). 
The major difference is that the square root is inside 
rather than outside the integral sign. 

The parametrised Jacobi action is of the form / = 
\TQ\f~K 1 \ where Q is the potential energy and K' is 
the kinetic energy, a function quadratic in the velocities. 
Therefor I consider a general parametrised action of the 
'local square root' form, i.e. 
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where a and /3 are again arbitrary constants. We can 
work out the momentum via 
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The momentum is homogeneous of degree zero in the ve- 
locity, which is the key property of a 'square root' action. 
This guarantees that the momenta are not independent. 
If I pick 



a 



1 

A' 



B 



ZAB - 2A 2 ' 
one can show just by direct substitution that 
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AM ab M cd P ac P bd - | [M ab P ab ] 2 



[detM] Q = 0. (24) 



This is a primary constraint, which is identical to the 
constraint which generated the Hamiltonian, Eq.©. It 
is easy to show that the Euler-Lagrange equations of 
the Lagrangian, Eci. (|20[) . exactly agree with Hamilton's 
equations of the Hamiltonian, Eq.JHJ). If we make the 
'good' choice of the constants A = B = 1, we get 
a = (3 = 1. This is the standard switch that occurs in go- 
ing from the ADM Lagrangian to the ADM Hamiltonian 
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where the kinetic energy goes from K ab K a b — (trK) 2 to 

A more sophisticated (but, I think, fairly dishonest) 
argument could be made as follows: If one has a solution 
to the Maxwell equations on a closed manifold without 
boundary, the total charge must vanish, from just the 
integral of Poisson's equation. Something similar might 
be expected to happen in gravity. Therefore the total 
gravitational energy should vanish which might give an 
integral constraint of the form J (K — Q') = 0, where K 
is the kinetic energy, quadratic in the momenta, and Q' 
is the negative of the potential energy (some function of 
the configuration variables). This means we can write 
K — Q' = divV where V is some vector. If this vector 
were a function of the configuration variables only, and 
it would have to be if we were not to have derivatives 
of the momenta, then we can absorb it into the Q via 
Q = Q' + divV. This promotes a global constraint to a 
local one of the desired form. 

Caveats: There are a number of (minor) warnings that 
should issued. As mentioned earlier, an undifferentiated 
term can be added to the potential which will just mimic 
a cosmological constant, there is no restriction on either 
magnitude or sign. I started with a three dimensional 
space and ended with four dimensional spacetime. This 
procedure works in any number of dimensions, all that 
changes is the value of the coefficient of the (M a i,P ab ) 2 
term. In n space dimensions it becomes — l/(n — 1). 
Finally, this approach does not pick up the signature of 
spacetime. At one point I used a scaling argument to set 
A = 1, C = 1. I could just as easily have set A = 1, C = 
— 1. This will generate Euclidean gravity. 

The approach advocated here can be considered as an 
enforced Yang-Mills construction. In Yang-Mills one is 
given a global symmetry, deliberately promotes it to a 
local symmetry and evaluates the consequences. Here 
the outcome is the same, a global symmetry becomes 



a local symmetry, but I had no choice in the matter. 
The requirement that the constrained Hamiltonian be 
selfconsistent forced the localization on me. 

This is certainly not the 'best' result of this kind. Fol- 
lowing I expect that if one added a scalar and/or 
vector(s) to the initial configuration space, the scalar 
wave equation, Maxwell theory, and Yang-Mills theory 
should emerge. One needs to look at the case where 
higher derivative terms appear in the potential. I have 
not done this, but based on what happened in Q, I would 
be surprised if anything new appeared. More interest- 
ingly, one could relax the condition that the momenta be 
quadratic and undifferentiated. It may be that something 
new emerges. Finally, one could drop the symmetry con- 
dition on the matrices. It may be that then one is led to 
something akin to the Cartan theory of gravity. 



* niall@ucc.ie !9ol 

[1] R. Arnowitt, S. Deser, and C. Misner in Gravitation: an 
Introduction to Current Research, ed. L. Witten, (Wiley, 
New York, 1962); C. Misner, K. Thorne, J. Wheeler, Grav- 
itation (Freeman, San Francisco, 1973). 

[2] S. Hojman, K. Kuchaf, C. Teitelboim, Ann. Phys (N.Y.) 
96, 88 (1976). 

[3] J. Barbour, B. Foster, N. O Murchadha, Class. Quant. 
Grav. 19, 3217 (2002). 

[4] N. 6 Murchadha, IJMP, 20, 2717 (2002). 

[5] P. Dirac, Lectures on Quantum Mechanics, (Yeshiva Uni- 
versity, New York, 1964). 

[6] C. Isham, Proc. Roy. Soc, A 351, 209 (1976); C. Teit- 
elboim, Phys. Rev. D 25, 3159 (1982); M. Pilati, Phys. 
Rev. D 26, 2645 (1982), ibid D 28, 729 (1983). 

[7] J. Barbour and N. 6 Murchadha, |gr-qc/99lfuT7l E. An- 
derson, J. Barbour, B. Foster, N. O Murchadha, to be 
published in Class. Quant. Grav. (2003), gr-qc/0211022 

[8] C. Lanczos, The Variational Principles of Mechanics, 
(Dover, New York, 1986). 



